PLBD: protein-ligand binding database of thermodynamic and kinetic intrinsic parameters.

We introduce a protein-ligand binding database (PLBD) that presents thermodynamic and kinetic data of reversible protein interactions with small molecule compounds. The manually curated binding data are linked to protein-ligand crystal structures, enabling structure-thermodynamics correlations to be determined. The database contains over 5500 binding datasets of 556 sulfonamide compound interactions with the 12 catalytically active human carbonic anhydrase isozymes defined by fluorescent thermal shift assay, isothermal titration calorimetry, inhibition of enzymatic activity and surface plasmon resonance. In the PLBD, the intrinsic thermodynamic parameters of interactions are provided, which account for the binding-linked protonation reactions. In addition to the protein-ligand binding affinities, the database provides calorimetrically measured binding enthalpies, providing additional mechanistic understanding. The PLBD can be applied to investigations of protein-ligand recognition and could be integrated into small molecule drug design. Database URL https://plbd.org/.

Structure and mechanism of secondary sulfonamide binding to carbonic anhydrases.

Zinc-containing metalloenzyme carbonic anhydrase (CA) binds primary sulfonamides with extremely high, up to picomolar, affinity by forming a coordination bond between the negatively charged amino group and the zinc ion and making hydrogen bonds and hydrophobic contacts with other parts of the inhibitor molecule. However, N-methyl-substituted, secondary or tertiary sulfonamides bind CA with much lower affinity. In search for an explanation for this diminished affinity, a series of secondary sulfonamides were synthesized and, together with analogous primary sulfonamides, the affinities for 12 recombinant catalytically active human CA isoforms were determined by the fluorescent thermal shift assay, stopped-flow assay of the inhibition of enzymatic activity and isothermal titration calorimetry. The binding profile of secondary sulfonamides as a function of pH showed the same U-shape dependence seen for primary sulfonamides. This dependence demonstrated that there were protein binding-linked protonation reactions that should be dissected for the estimation of the intrinsic binding constants to perform structure-thermodynamics analysis. X-ray crystallographic structures of secondary sulfonamides and computational modeling dissected the atomic contributions to the binding energetics. Secondary sulfonamides bind to carbonic anhydrases via coordination bond between the negatively charged nitrogen of alkylated amino group and Zn(II) in the active site of CA. The binding reaction is linked to deprotonation of the amino group and protonation of the Zn(II)-bound hydroxide. To perform the structure-thermodynamics analysis, contributions of these linked reactions must be subtracted to determine the intrinsic energetics. In this aspect, the secondary sulfonamides are similar to primary sulfonamides as CA inhibitors.

Standard operating procedure for fluorescent thermal shift assay (FTSA) for determination of protein-ligand binding and protein stability.

A standard operating procedure for a fluorescence-based thermal shift assay (FTSA) is provided describing its typical applications, advantages and limitations. FTSA is a simple, robust, universal and quick assay to determine protein-ligand binding affinities and protein stabilities in the presence of various excipients and solution conditions. Therefore, the assay is very useful for the straightforward characterization of new recombinantly produced proteins. The assay has a wide dynamic range enabling simultaneous determination of affinities in the milimolar to picomolar range. The assay could be used for essentially any protein that is sufficiently soluble and stable in the studied aqueous solution. Here we provide examples and typical experimental protocols for both affinity and stability determinations.

Uncertainty in protein-ligand binding constants: asymmetric confidence intervals versus standard errors.

Equilibrium binding constants (Kb) between chemical compounds and target proteins or between interacting proteins provide a quantitative understanding of biological interaction mechanisms. Reported uncertainties of measured experimental parameters are critical for decision-making in many scientific areas, e.g., in lead compound discovery processes and in comparing computational predictions with experimental results. Uncertainties in measured Kb values are commonly represented by a symmetric normal distribution, often quoted in terms of the experimental value plus-minus the standard deviation. However, in general, the distributions of measured Kb (and equivalent Kd) values and the corresponding free energy change ΔGb are all asymmetric to varying degree. Here, using a simulation approach, we illustrate the effect of asymmetric Kb distributions within the realm of isothermal titration calorimetry (ITC) experiments. Further we illustrate the known, but perhaps not widely appreciated, fact that when distributions of any of Kb, Kd and ΔGb are transformed into each other, their degree of asymmetry is changed. Consequently, we recommend that a more accurate way of expressing the uncertainties of Kb, Kd, and ΔGb values is to consistently report 95% confidence intervals, in line with other authors' suggestions. The ways to obtain such error ranges are discussed in detail and exemplified for a binding reaction obtained by ITC.

Halogenated and di-substituted benzenesulfonamides as selective inhibitors of carbonic anhydrase isoforms.

By applying an approach of a "ring with two tails", a series of novel inhibitors possessing high-affinity and significant selectivity towards selected carbonic anhydrase (CA) isoforms has been designed. The "ring" consists of 2-chloro/bromo-benzenesulfonamide, where the sulfonamide group is as an anchor coordinating the Zn(II) in the active site of CAs, and halogen atom orients the ring affecting the affinity and selectivity. First "tail" is a substituent containing carbonyl, carboxyl, hydroxyl, ether groups or hydrophilic amide linkage. The second "tail" contains aryl- or alkyl-substituents attached through a sulfanyl or sulfonyl group. Both "tails" are connected to the benzene ring and play a crucial role in selectivity. Varying the substituents, we designed compounds selective for CA VII, CA IX, CA XII, or CA XIV. Since due to binding-linked protonation reactions the binding-ready fractions of the compound and protein are much lower than one, the "intrinsic" affinities were calculated that should be used to study correlations between crystal structures and the thermodynamics of binding for rational drug design. The "intrinsic" affinities together with the intrinsic enthalpies and entropies of binding together with co-crystal structures were used demonstrate structural factors determining major contributions for compound affinity and selectivity.

Intrinsic Thermodynamics of Protein-Ligand Binding by Isothermal Titration Calorimetry as Aid to Drug Design.

Isothermal titration calorimetry (ITC) is one of the main techniques to determine specific interactions between molecules dissolved in aqueous solution. This technique is commonly used in drug development programs when low-molecular-weight molecules are sought that bind tightly and specifically to a protein (disease target) molecule. The method allows a complete thermodynamic characterization of an interaction, i.e., ITC enables direct determination of the model-independent observed interaction change in enthalpy (ΔH) and a model-dependent observed interaction affinity (change in Gibbs free energy, ΔG) in a single experiment. The product of temperature and change in entropy (TΔS) can be obtained by the subtraction of ΔG from ΔH, and the change in heat capacity (ΔCp) can be determined as a slope of the temperature dependence of the binding ΔH.Despite the apparent value of ITC in characterization of interactions, it is often forgotten that many protein-ligand binding reactions are linked to protonation-deprotonation reactions or various conformational changes. In such cases, it is important to determine the linked-reaction contributions and obtain the intrinsic values of the changes in Gibbs energy (affinity), enthalpy, and entropy. These energy values can then be used in various SAR-type structure-thermodynamics and combined with structure-kinetics correlations in drug design, when searching for small molecules that would bind the protein target molecule. This manuscript provides a detailed protocol on how to determine the intrinsic values of protein-ligand binding thermodynamics by ITC.

Repeatability, precision, and accuracy of the enthalpies and Gibbs energies of a protein-ligand binding reaction measured by isothermal titration calorimetry.

In rational drug design, it is important to determine accurately and with high precision the binding constant (the affinity or the change in Gibbs energy, ∆G), the change in enthalpy (ΔH), and the entropy change upon small molecule drug binding to a disease-related target protein. These thermodynamic parameters of the protein-ligand association reaction are usually determined by isothermal titration calorimetry (ITC). Here, the repeatability, precision, and accuracy of the measurement of the affinity and the change in enthalpy upon acetazolamide (AZM) interaction with human carbonic anhydrase II (CA II) are discussed based on the measurements using several ITC instruments. The AZM-CA II reaction was performed at decreasing protein-ligand concentrations until the determination of ∆G and ΔH was not possible, indicating a lower limit for accuracy. To obtain the confidence intervals (CI) of the ∆G and ΔH of AZM binding to CA II, the binding reaction was repeated numerous times at the optimal concentration of 10 µM and 25 °C temperature. The CI (at a confidence level α = 0.95) for ΔH = - 51.2 ± 1.0 kJ/mol and ∆G = - 45.4 ± 0.5 kJ/mol was determined by averaging the results of multiple repeats.

Pyrrolidinone-bearing methylated and halogenated benzenesulfonamides as inhibitors of carbonic anhydrases.

Two series of benzenesulfonamides bearing methyl groups at ortho/ortho or meta/ortho positions and a pyrrolidinone moiety at para position were synthesized and tested as inhibitors of the twelve catalytically active human carbonic anhydrase (CA) isoforms. Observed binding affinities were determined by fluorescent thermal shift assay and intrinsic binding affinities representing the binding of benzenesulfonamide anion to the Zn(II)-bound water form of CA were calculated. Introduction of dimethyl groups into benzenesulfonamide ring decreased the binding affinity to almost all CA isoforms, but gained in selectivity towards one CA isoform. A chloro group at the meta position of 2,6-dimethylbenzenesulfonamide derivatives did not influence the binding to CA I, but it increased the affinity to all other CAs, especially, CA VII and CA XIII (up to 500 fold). The compounds may be used for further development of CA inhibitors with higher selectivity to particular CA isoforms.

Thiazole-substituted benzenesulfonamides as inhibitors of 12 human carbonic anhydrases.

Four series of para or meta - substituted thiazolylbenzenesulfonamides bearing Cl substituents were designed, synthesized, and evaluated as inhibitors of all 12 catalytically active recombinant human carbonic anhydrase (CA) isoforms. Observed affinities were determined by the fluorescent thermal shift assay and the intrinsic affinities were calculated based on the fractions of binding-ready deprotonated sulfonamide and CA bearing protonated hydroxide bound to the catalytic Zn(II) in the active site. Several compounds exhibited selectivity towards CA IX, an anticancer target. Intrinsic affinities reached 30 pM, while the observed affinities - 70 nM. The structure-intrinsic affinity relationship map of the compounds showed the energetic contributions of the thiazole ring and its substituents.

Benzimidazole design, synthesis, and docking to build selective carbonic anhydrase VA inhibitors.

The similarity of human carbonic anhydrase (CA) active sites makes it difficult to design selective inhibitors for one or several CA isoforms that are drug targets. Here we synthesize a series of compounds that are based on 5-[2-(benzimidazol-1-yl)acetyl]-2-chloro-benzenesulfonamide (1a) which demonstrated picomolar binding affinity and significant selectivity for CA isoform five A (VA), and explain the structural influence of inhibitor functional groups to the binding affinity and selectivity. A series of chloro-substituted benzenesulfonamides bearing a heterocyclic tail, together with molecular docking, was used to build inhibitors that explore substituent influence on the binding affinity to the CA VA isoform.

N-Sulfamoylphenyl- and N-sulfamoylphenyl-N-thiazolyl-ß-alanines and their derivatives as inhibitors of human carbonic anhydrases.

A series of N-substituted and N,N-disubstituted ß-amino acids and their derivatives bearing benzenesulfonamide moiety were designed and synthesized in search of compounds that would be high-affinity and selective inhibitors of human carbonic anhydrases (CA). There are 12 catalytically active human CA isoforms, the cytosolic CA I, CA II, CA III, CA VII, and CA XIII, secreted CA VI, the mitochondrial CA VA and CA VB, membrane-associated CA IV, and transmembrane CA IX, CA XII, and CA XIV. The di-bromo meta-substituted compounds exhibited low nanomolar dissociation constants and over 10-fold selectivity for mitochondrial isozyme CA VB, implicated in diseases of the central nervous system and obesity. These compounds can be used for further development as inhibitors of significant binding affinity and selectivity towards CA VB isozyme.